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ABSTRACT 

We study the capability of the James Webb Space Telescope (JWST) to detect 
Supermassive Dark Stars (SMDS). If the first stars are powered by dark matter heating 
in triaxial dark matter haloes, they may grow to be very large > lO^M© and very 
bright > 10^ Lq. These SMDSs would be visible in deep imaging with JWST and 
even Hubble Space Telescope (HST). Indeed the object detected at z '-^ 10 in the 
recent HST ultra deep field image with the Wide Field Camera 3 (WFC3) may be a 
SMDS. We use sensitivity limits from previous HST surveys to place bounds on the 
numbers of SMDSs that may be detected in future JWST imaging surveys. We showed 
that SMDS in the mass range 10^ — lO'^M© are bright enough to be detected in all 
the wavelength bands of the NIRCam on JWST (but not in the less sensitive MIRI 
camera at higher wavelengths). If SMDSs exist at z ~10, 12, and 14, they will be 
detectable as J-band, H-band, or K-band dropouts, respectively. With a total survey 
area of 150 arcmin^ (assuming a multi-year deep parallel survey with JWST), we 
find that typically the number of IO^AIq SMDSs found as H or K-band dropouts is 
^ 10^ fsMDS 1 where the fraction of early DM haloes hosting DS is likely to be small, 
fsMDS << 1- If the SDMS survive down to z=10 where HST bounds apply, then the 
observable number of SMDSs as H or K-band dropouts with JWST is ^ 1 — 30. While 
individual SMDS are bright enough to be detected by JWST, standard PopIII stars 
(without dark matter annihilation) are not, and would only be detected in first galaxies 
with total stellar masses of 10^ — IO^Mq. Differentiating first galaxies at z>10 from 
SMDSs would be possible with spectroscopy: the SMDS (which are too cool produce 
significant nebular emission) will have only absorption lines while the galaxies are likely 
to produce emission lines as well. Of particular interest would be the Hell emission 
lines at A 1.6/i as well as Ha lines which would be signatures of early galaxies rather 
than SMDSs. The detection of SMDSs with JWST would not only provide alternative 
evidence for WIMPs but would also also provide a possible pathway for the formation 
of massive (10^ — IO^Mq) seeds for the formation of supermassive black holes that 
power QSOs at z = 6. 

Key words: dark matter - first stars - stars: Population HI - stars: pre-main- 
sequence-galaxies: high-redshift 



1 INTRODUCTION 

The first stars are thought to have formed at z = 10 — 50 
when the universe was about 200 million years old in ~ 
10^ Mq (mini) halos consisting of 85% DM and 15% baryons 
in the form of H and He from big bang nucleosynthesis. Their 
formation marks the end of the "dark ages" of the Uni- 
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verse. For reviews o f the standard picture o f the formation 
of the first stars see iBarkana fc Loebl ll200ll); lYoshida et al.l 



(2003'):'Bromm fc LarsonI (|2004l '): lRipamonti fc Abell (120051 ): 
.Brom m ct al. (20091). 

ISpolv ar et al. I (|2008t ) first showed that dark matter 
heating may drastically alter the picture of formation for 
these first stars. The canonical example of particle DM is 
Weakly Interacting Massive Particles (WIMPs). In many 
theories WIMPs are their own antiparticles and annihilate 
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with themselves wherever the DM density is high. In fact, 
this annihilation process is exactly what is responsible in the 
early Universe for leaving behind the correct relic WIMP 
abundance today to solve the dark matter problem, 24% of 
the energy density of the universe. The same annihilation 
process would then take place also in the collapsing proto- 
stellar clouds at the centers of minihalos. At suitably high 
baryonic density in these clouds, the annihilation products 
get stuck inside the cloud and prevent it from undergoing 
further collapse. The annihilation products thermalize with 
the baryons and provide a very powerful heat source. In- 
deed, the object becomes a "dark star", which, despite its 
name, shines very bright. The DM - while only a negligible 
fraction of the star's mass - provides the key power source 
for the star through DM heating. Note that the term 'Dark' 
refers to the power source, not the content of the star. These 
first Dark Stars are stars made primarily of hydrogen and 
helium with a smattering of dark matter (<1% of the mass 
consists of DM); yet they shine due to DM heating. 

Recently there has been much excitement in the 
dark matter community about hints of WIMP detec- 
tions in a number of e xperiments: excess positrons in 
the P AMELA sate l lite jAdriani et all l2009l : lAbdo et all 
I2OIOI : lAdriani etakl I2OI0I ') may be due to DM annihila- 
tion (though alternative astrophysical explanati ons are more 
likely) . Excess 7-rays in th e FERMI satellite (Abdo et al.' 
2009bl lal: iDobler et all I2OI0I : iThe Fermi LAT Collaboration 



20 111 ) m ay be due to DM annihilation; and an nual mod- 



ulation (jDrukier et al.l 1 19861 : iFreese et al. Il988l ) in direct 
detection e xperiments DAMA ( Bernabei et al.l I2OI0I ) and 
COGENT llAalseth et aLll201ll ). The CRESST experiment 
ijAngloher et al.l 2011 ) also has unexplained events. 



The WIMP annihilation rate is n^{av) where is 
WIMP density and we take the standard annihilation cross 
section 



(av) = 3 X 10" 



s 3 / 
cm /s, 



(1) 



and WIMP masses in the range 1 GeV-10 TeV. WIMP an- 
nihilation produces energy at a rate per unit volume 



Qdm = n^{av)my, = {av)p^/m^ 



(2) 



where is the WIMP number density, is the WIMP 
mass, and is the WIMP energy density. The annihila- 
tion products typically are electrons, photons, and neutri- 
nos. The neutrinos escape the star, while the other anni- 
hilation products are trapped in the dark star, thermalize 
with the star, and heat it up. The luminosity from the DM 
heating is 



[dV 



(3) 



where fg is the fraction of the annihilation energy deposited 
in the star (not lost to neutrinos) and dV is the volume 
element. We take /q = 2/3 as is typical for WIMPs. 

Dark stars are born with masses ~ IMq . They are giant 
puffy (~ 10 AU), c ool (surface temper atures < 10, OOOK), 
yet bright objects (jFreese et al.l l2008al ). They reside in a 
large reservoir (~ 10^ M©) of baryons, i.e., ~ 15% of the 
total halo mass. These baryons can start to accrete onto 
the dark stars. Dark stars can continue to grow in mass 
as long as there is a supply of DM fuel. We consider two 
different mechanisms that can continually provide the req- 



uisite dark matter fuel, allowing them to become supermas- 
sive dark stars (SMDS) of mass Mds > IO^A/q. : 

1) Extended Adiabatic Contraction (AC): the central 
dark matter density is enhanced due to an increase in the 
depth of the gravitation potential well due to the infall of 
baryons. We treat this gravitational effect via the Blumen- 
thal method for adiabatic c ontraction. While th is approach 
is simple to imp l ement , we Freese et a l. (2009) and others 
iNataraian et al.l (|2009l ):ll occo et al.l (|2008l have previously 
shown that it provides dark matter densities accurate to 
within a factor of two, which is perfectly adequate for these 
studies. In the central cusps of triaxial DM halos DM par- 
ticles follow a variety of cen trophilic orbits (box orbits and 
chaotic orbits FValluri et al.lfiol O') whose population is con- 
tinuously replenished, allowing DM annihilation to continue 
much longer than in spherical DM halos. The period of ex- 
tended AC can thus last for a very long time (h undreds 
of millions of years or more). iFreese et al] (|2010al ) showed 
that this replenishment of the DM in the central cusp could 
be used to followed the growth of dark stars from their in- 
ception at IMq, till they become supermassive dark stars 
(SMDS) of mass Mds > 10^ Mq. 

2) Capture: As a second mechanism for dark matter re- 
fueling, we take the star to be initially powered by the 
DM from adiabatic contraction (AC), but assume the AC 
phase is short 300, 000 years; once this DM runs out, 
the star shrinks, its density increases, and subsequently 
the DM is replenished i nside the star by capture of DM 
from the surroundings (jFreese et al.l l2008bl : lloccd I2OO8I : 
ISivertsson fc Gondold[201ll ) as it scatters elastically off of 
nuclei in the star. In this case, the additional particle physics 
ingredient of WIMP scattering is required. This elastic scat- 
tering is the same mechanism that direct detection exper- 
iments (e.g. CDMS, XENON, LUX, DAMA, COGENT, 
COUPP, CRESST) are using in their hunt for WIMPs. 

Supermassive dark stars can result from either of these 
mechanisms for DM refueling inside the star. lUmeda et al.l 
(2009) considered a different scenario which also results in 
SMDSs. In all of these cases SMDSs can live for a very long 
time, tens to hundreds of million years, or possibly longer 
(even to today). We find that SMDS of mass Mds > W'^Mq 
SDMSs are very bright > 3 x 10^ Lq which makes them 
potentially observable by the James Webb Space Telescope 
(JWST). ' 

The key ingredient that allows dark stars to grow so 
much larger than ordinary fusion powered Population III 
stars is the fact that dark stars are so much cooler. Ordinary 
Pop III stars have much larger surface temperatures in ex- 
cess of 50, OOOK. They produce ionizing photons that provide 
a va riety of feedback mech anisms that cut off further accre- 
tion. [McKeeA^^ (120081 ) have estimated that the resultant 
Pop III stellar masses are ~ I4OM0 . The issue of the initial 
mass function for Po p III stars is far from being solved. R e- 
cent simulations Csee iClark et al.ll201ll : iGreff et al.ll2011al lbl) 
indicate that the typical mass of such objects is much lower 
that previously thought. Dark stars are very different from 
fusion-powered stars, and their cooler surface temperatures 
allow continued accretion of baryons all the way up to enor- 
mous stellar masses, Mds > 10^5 M©. 

In this paper we discuss detectability of these objects 
in the upcoming JWST. In future work we will investigate 
how well other observations with Herschel, SPITZER, GMT, 
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TMT and other instruments can detect or place bounds on 
Dark Stars. We restrict our discussion only to SMDS of 
mass 10'^, lO^M© (we show that SM DS of lO '^ and lower are 
hard to detect). Previously Zackriss on et al.l (|2010a ) studied 
dark stars of even lower masses, and concluded that even 
<10^Mq DS could be detected as individual objects with 
JWST if their fluxes were magnified by gravitational lens- 
ing by a well-placed foreground cluster. Since supermassive 
dark stars are larger and brighter, they are easier to de- 
tect. A preliminary study of detectability wit h JWST and 
HST o f sup ermassive Dark Star s was made in iFreese et al.l 
(l2010al ) an d iFreese et al](j2010bl ). Freese et al. approximated 
the spectrum of the SD as a pure blackbody determined by 
its temperature and radius and used it to show that indi- 
vidual SMDS would be detectable with JWST and HST. 
In this paper we improve our estimate by using spectra 
from the TLUSTY model stellar atmosphe res code for zero- 
met a llicity atmospheres from the work of IZackrisson et"al] 

HoTobi). 

SMDS formed via Extended AC are easier to detect 
than those formed with capture. Those formed "with cap- 
ture" are somewhat hotter (by a factor of few) and have 
radii smaller by a factor of 5-10 for the same stellar mass. 
Because they are hotter, their peak wavelength moves out 
of the most sensitive ranges for HST and JWST, and their 
fluxes in the detectors are lower. 

Once the SMDS run out of DM fuel, they contract and 
heat up till the core reaches 10*K and fusion begins. Due 
to their extremely large masses the fusion-powered phase 
is short and the SMSDs collapse to from massive black 
holes of mass 10* - 10*^ Mq. a gain, this prediction is dif- 
ferent from standard Pop HI stars, many of which ex plode 
as pair-instabilty supernovae (|Heger fc Wooslevll2002l ) with 
predicted even/odd element abundance ratios that are not 
(yet) observed in nature. These massive black holes rem- 
nants could provide the moderately massive "seeds" for the 
formation of nuclear supermas sive black holes accoun ting for 
the existence of lO" M© BHs (|Haiman fc LoeblbOOll ) which 
are the central engines of the most distant {z > 6) quasars in 
the S loan Digital Sky Survey (SDSS) (Fa n et al. 2001. 2004,, 
l2006l ). Indeed direct collapse of very metal-po or, low-angular 
momentum gas via dynamical instabilities (iLoeb fc Rasid 
1 19941 : iBegelman et al.ll2006l : iLodato fc Nataraianlbood ) has 
been proposed as a way to form massive "seed" black holes 
of 10"* - IQ^Mq at redshifts of 10-15. These massive seed 
formation scenarios however, are difficult to confirm obser- 
vationally since the BHs form in compact, low luminosity 
cold gas disks and the BH formation is accompanied by a 
sudden burst of with a luminosity of W^Lq. In contrast if the 
"seeds" form from SMDSs, they may well shine for 10® — 10^ 
years prior to their collapse to a BH, enabling them to be 
detected by JWST. 

SMDS could also make plausible precursors of Interme- 
diate Mass Black Holes; and account for the BHs inferred 
by ext ragalactic radio exc ess seen by the ARCADE exper- 
iment jSeiffert et al.ll2009l ). In addition, the BH remnants 
from DS could play a role in high-redshift gamma ray bursts 
thought to take place due to accretion onto early black holes 
iNaravan et al.ll200ll ') 

The possibility that DM annihilation might have 
effects on today 's stars was initiall y considered in the '80s 
and early '90s (|Krauss et all Il985l : iPress fc Spergell 1 19851 : 



iBouauet fc Salatilll989l : ISalati fc Silklll989l) and has recently 
been studi ed in intere stin g pap ers b y iMoskalenko fc Wail 



(2007 
(200a), 



IScott et al 



Scott et al. 



(120071'): i Bertone 



Scotj 120101 ) 



Fairbaini 
Casanellas fc Lopes ( 2009 



iHooper et al.l (|201C 

Several authors have explored the re percussions of 
DM h ea ting in the fi rst sta rs i ncluding: ISpolvar et al.] 



iTaoso et al 



(l2008l): iFreese et all ll2008bllal): ISpolvar et al 



Yoon et al.l 



iRipamonti et al.l 



locco et al.l 
I Schleicher et al.l 



(•2010); 'Ripamon ti et al 
(2011); Casanellas fc Lopea 




Scotd 



'Gondolo et al. 
[Sivcrtsson fc Gondolo 

Hira no et all (|201ll ): lloccol (|2010l ): lllie et al] (|201ll '):" 
{201^. 

The effects of DS (and those of the re sultant main se- 
quence stars) on reionization was s tudied b y'Schleiche r et al] 
(2008, 2009) and more recently by Scott et al. (2011) as dis- 
cussed below. 

In this pa p er we follow the approach taken by 
IZackrisson et al.l (|2010bl . [20TTbl 'l. Similar to their work, we 
use SMDS spectra from the TLUSTY code; compute the for- 
mation rate of DSs by counting DM haloes in N-body simu- 
lations; and use HST data to bound the numbers of SMDS 
that survive to z=10 and therefore the numbers that may 
be seen with JWST. Their study focused on IO^'Mq SMDS 
while we consider lower mass ones as well. We go beyond 
their work by studying SMDS as H and K-band dropouts 
with JWST, where JWST can really improve upon all pre- 
vious data sets. 

The paper is organized as follows. In Section II we de- 
scribe the spectra of SMDS obtained by using the TLUSTY 
code. In Section HI we compute formation rate of DSs, by 
counting DM halos in a N-body simulation of structure for- 
mation at z > 10 carried out with the Cube P^M code 
(|lliev et al.ll2010l ) and assuming that some fraction fsMDS of 
these early halos will host DS. In Section IV, we examine the 
detectability of SMDS in Hubble Space Telescope (HST). In 
fact HST has seen objects out to 2: ~ 10, and it is interesting 
to speculate that HST could already have seen SMDS if they 
survive to redshift z=10. With current imaging data it is im- 
possible to differentiate between an early galaxy composed of 
PopIII stars from a SMDS. However, the the fact that HST 
has only seen one object at this high redshift can be used 
to set bounds (IZackrisson et alll2010al ) the numbers of dark 
stars at 2 ~ 10. In Section V we show that dark stars may 
be detected in a variety of JWST filters, and in particular 
may show up as J-band, H-band, or K-band dropouts; such 
a detection would then gives an indication of their redshift. 
In Section VI, we compare early galaxies at high redshifts 
(consisting of PopIII stars with different IMFs) with SMDS, 
which will look very similar with JWST, and start a discus- 
sion of ways to differentiate between them. In Section VII 
we conclude and summarize the results of our study. 



2 DARK STAR SPECTRA 

In this section we present spec tra of SMDS ob tained with the 
pubhcly available TLUSTY (|Hubenvlll988l') svnth e tic ste l- 
lar atmospheres code. As discussed in IFreese et al.l (|2010al ). 
SMDS formed via captured DM are much hotter than SMDS 
formed via extended AC. Also, stars formed via capture 
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undergo a Kelvin-Helmholz contraction phase prior to DM 
capture, hence their radii are 5-10 times smaller than those 
SMDS ol the same mass formed via the extended AC mech- 
anism. Since Dark Stars are composed of primordial hy- 
drogen and helium, no other elements are assumed to be 
present in the atmosphere and, hence all the observed spec- 
tral lines are those of H and He. However the differences in 
the temperature and radii of SMDS formed via these two 
mechanisms are responsible for the differences in the spec- 
tra in the two panels of Figure [T] The left panel shows 
the spectrum for a l(f Mq dark star with surface tempera- 
ture Te// — 1.9 X 10*K which grew via extended AC. The 
Lyman edge is seen at roughly 0.1 micron^. Similarly, Fig- 
ure [T] (right ) illustrates the spectrum for a IO'^Mq and 
Te// = 5.1 X lO^'K DS which grew via captured DM. The 
most prominent differences from the left panel are a shift of 
the peak in the spectrum to lower wavelengths and a steeper 
UV continuum slope P {fx oc A^). Despite the fact that the 
SMDS formed via capture is hotter, its significantly smaller 
radius makes it harder to detect in the near infrared at red- 
shifts of ~ 10 and above. 

There are significant differences in the spectra the two 
cases. In the left panel (extended AC), the lower surface 
temperature (~ 2 x lO'^K) implies a significant fraction of 
neutral H and He remain in the stellar atmosphere result- 
ing in strong absorption lines at wavelengths correspond- 
ing to the Lyman series (0.1216 )i.m-0.0912 \im). At shorter 
wavelengths we notice another break in the spectrum due 
to neutral helium (Hel) absorption (~0.05 nm-0.06 vim). In 
the right panel ("with capture"), the higher surface temper- 
ature (Te// ~ 5 X 10* K) implies that H is ionized hence the 
Lyman absorption lines are weaker. The break in the spec- 
trum in Figure [T] right panel corresponds to absorption by 
singly ionized helium (Hell) at wavelengths ranging between 
0.023 [im and 0.030 \im. In the left panel Hel lines appear at 
wavelengths ~ [0.3 (j.m-0.45 vim]. Hell lines at wavelengths 
~0.46 vim, and more Hel lines at ~ [0.47 vim-0.7 vim]. The 
same lines, with somewhat weaker strength, are seen in the 
right panel. In both cases we note a sequence of absorption 
lines between ~[0.8 vim- 1.0 vim], which correspond to Hel 
absorption. 



3 DARK STAR FORMATION RATE 

The first Dark Stars can form in the early Universe inside 
minihalos of ~ lO^M©, where protostellar clouds collapse 
via molecular hydrogen cooling until the DM heating sets 
in. Later in 10* Mq halos, where clouds collapse via atomic 
line radiative cooling, larger DS can form. To compute the 
detection rate of SMDS with JWST we need to know the 
formation rate of 10^ — 10* dark matter halos. If we as- 
sume that a fraction fsMDS of these halos contain Dark Stars 
we can use this to compute the formation rate of DSs. We 

^ Compared to a blackbody of the same temperature, photons 
below the Lyman edge have typically been shifted to higher wave- 
lengths (lower energy) by absorption and rescattering. However, 
the excess seen at wavelengths just below the Lyman edge is due 
to photons coming from deeper inside the star (the photosphere 
is at roughly an optical depth ~ 1, and at this wavelength there 
is very little absorption). 



will attempt to set constraints on this fraction by using the 
fact that a single z = 10 object was observed in recent HST 
ultra deep field obser vations with the Wide Field Camera 3 
(jBouwens et al.ll201ll . h ereafter HUDF09 ). 

A similar study bv IZackrisson et all l|2010bl . l2011bl ) for 
the case of lO^M© SMDS concluded that the prior null de- 
tection of z = 10 obj ects in first year HUDF09 observation 
(jBouwens et al.ll20lrf ). was sufficient to rule out the detec- 
tion of 10'' M0 SMDS with JWST. However these authors 
did not consider the effect of the time it takes the SMDS to 
grow when computing the formation rates for DM halos that 
could host such objects. This effect is transparent in Table[T] 
in the differences between what we labeled as Zstart (the red- 
shift that should be used when computing the formation rate 
of DM halos) and Zform (the redshift when the DS reaches 
its final mass). Consideration of a finite time required for 
the SMDS to grow (following the formation of its host DM 
halo) significantly lowers the bounds predicted from HST, 
since to be visible at z = 10 the more massive DM halos has 
to have formed at a higher redshifts, where they are rarer. 
In addition we consider the case of the 10^ Mq SMDS, since 
these objects are likely to be more numerous, are detectable 
with JWST, and are also subject to bounds from existing 
HST observations. 

We use N-Bod y simulatio n s of s tructure formation at 
high redshifts from llhev et all HoiO) carried out with the 
CubeP^M N-Body code, dev eloped from the particle-mesh 
PM-FAST l|Merz et al.ll2005l ). This high resolution simula- 
tion considers a comoving volume of 6.3/i"^ Mpc with 1728^ 
particles of mass 5.19x IO^Mq, hence is able to resolve halos 
of mass > 5x IO^Mq. We compute the formation rate {dn/dt 
as a function of redshift per comoving Mpc'^ per year), of 
minihalos with masses within different mass ranges. Fig- 
ure [2] shows the formation rate of halos in two mass ranges 
that spans a factor of two in mass (lO'^ — 2 x IO'^Mq and 
10* — 2 X 10* Af©) while Figure |3] shows the formation rate of 
halos in two mass ranges that span a factor of five in mass 
(lO'^ - 5 X lO^M© and 10* - 5 x 10* Mq). 

We computed the formation rate of DM minihalos us- 
ing two different sets of bin widths, to show that the re- 
sults are relatively insensitive to this issue. As our canonical 
case, we computed the formation rate dn/dt of minihalos per 
Mpc~^yr~^ formed in a bin whose width is a factor of two 
in mass (Figure [2]). dn/dt as a function of redshift is shown 
for halos in the mass range (1 — 2) x IO'^Mq (left panel) and 
for halos in the mass range (1 — 2) x 10* Mq (right panel). We 
assume that since the baryonic fraction initially in the halo 
is roughly 15%, we assume that a DS forming in a halo of a 
given mass can attain at most 10-15% of the mass of its host 
halo. Following Freese et al. 10 we assume that the DS can 
grow with an accretion rate of ~ IMq to the point where it 
consumes a significant fraction of the baryons in the halo. 
In other words, we assume that a 10'' Mq SMDS wiU form 
in a (1 — 2) X 1O*M0 minihalo. While this is an unlikely sce- 
nario, which involves most of the baryons in the halo being 
accreted into a single central object, we will see that even 
with this assumption, detection rates of SMDS with JWST 
are fairly small. The formation rates in Scenario I are plotted 
in Figure [2] 

As a check, we also broadened the range of DM halo 
masses in which DS from by allowing halo masses to span 
a factor of five in mass. Figure |3] (left panel) shows the 
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Figure 1. Expected spectral energy distribution (SEDs) of IO^Mq supermassive dark stars. Left panel: DS with a surface temperature 
of 1.9 X 10*K and formed via the extended adiabatic contraction (AC) only mechanism. Right panel: with a surface temperature of 
5.1 X lO** and formed "with capture". 
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Figure 2. Left: the formation rate of 1 — 2 X IO'^Mq minihalos per comoving Mpc^ and year. These halos are potential hosts for the 
10^ Mq SMDS. Right: formation rate for 1 - 2 X 10® Mq minihalos in which a lO^'M© SMDS can form. The black lines correspond to the 
formation rate computed directly from the N-body simulation and the smoother red lines (obtained by computing a running average) 
simply improve visibility of the general trend. 



formation rate dn/dt as a function of redshift for halos in the 
mass range (1 — 5) x IO'^AIq and the right panel indicates the 
formation rate of halos in the mass range (1-5)x10**Mq. In 
this scenario the SMDS is 10 — 50 times smaller than its host 
halo, and is more realistic since in this case all the baryons 
in the halo are not accreted by the DS. 

A comparison of Figures [5] and [3] shows that the forma- 
tion rate of host halos does not vary significantly between 
the two scenarios (at most by a factor of 3) . Henceforth, in 
the remainder of this study, we will always take the halo 
mass range to span a factor of two in mass. 

We define Zstart to be the (approximate) formation red- 
shift of minihalos capable of hosting DS, allowing for an 



uncertainty of a unit redshift interval; i.e. the minimum 
redshift of minihalo formation is Zmin ~ Zstart — 1/2 while 
the maximum redshift is Zmax = Zstart + 1/2. We make a 
distinction between Zstart the redshift of formation of the 
DM halo capable of hosting a DS (initial ~ IMq mini 
dark stars come into existence very soon after this redshift), 
and Zform the redshift of formation of the SMDS. Between 
Zstart and Zform the DS grows by accreting baryons at a 
rate of 10~^ — 10~^MQ/yr growing over this period to a su- 
permassive size of - 10^ - 10'' Mq. This difference between 
Zstart and Zform is crucial to accounting for the differences 
between the results presented in this paper and previous 
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Figure 3. As in Figure [2] but with the larger mass bin width (see scenario II in text). In the left panel we plot the formation rate 
of minihalos with a mass in the 1 — 5 X 10^ Mq range, where a DS of 10® A/0 could form. The panel on the right is for halos in the 
1 - 5 X 10* M0 range, where a DS of IO'^Mq could form. 



work (jZackrisson et all l2010bh where this additional time 
required to grow supermassive was not allowed for. 

The formation rate of minihalos per unit redshift and 
arcmin^ is then given by 

-7-^ = "77 (Vc(zmin) - Vc(zmax)) ^ At{min; mux) , (4) 
azdO^ at iir 

where, Vc denotes the comoving volume at a given red- 
shift, C is the conversion factor between arcmin^ and stera- 
dians, and At{min; max) is the cosmic time interval between 

Zmin and Zniax- 

At{Tnin;max) — tn / -■ 

(5) 

TheN-body simulations from with the halo formation rates 
are computed as well as other calculations assume a stan- 
dard ACDM Universe in which Q,rn = 0.27 is the cosmic 
matter density and f^A = 0.73 is the cosmic dark energy den- 
sity or cosmological constan t with parameters from WMAP5 
data lKomatsu et all l|2009l '). 

We consider three possible redshifts Zform by which the 
DS has accreted enough baryons to become supermassive: 

• Case A: Zform =10 

• Case B: Zform =12 

• Case C: zform =15 

In principle the accretion rate and the final mass of the 
SMDS in these three cases will imply three different values of 
Zstart at which the relevant minihalos formed. To simplify the 
situation we assume a fixed accretion rate of W^AIq /yr to 
determine the values for Zatart (Table [T] column 5) and using 
corresponding dn/dt values from Figure [2] (Table [1] column 
6) at z =Zstart , we evaluate using Equation© (Table 

[1] column 7). These values of ^^^^ will be used in sections 
to follow. 



4 SUPERMASSIVE DARK STARS WITH 
HUBBLE SPACE TELESCOPE 

In this section we examine the observability of dark stars 
with existing Hubble Space Telescope (HST) surveys, spec- 
ulating that HST may already have seen such objects, if 
they survive to redshift z = 10. W e will adopt the stan dard 
"dropout technique" pioneered by ISteidel et al.l (|l996l ) and 
applied recently to J and H band observati ons of the Hub- 
ble Ultra Deep F ield (HUDF09) by (Bou wens et all 1201 ll : 
lOesch et al.ll201ll ) to detect a candidate galaxy at z = 10 
as a "J-band dropout." This photometric redshift determi- 
nation method requires a 5-sigma detection of an object in 
one band but a non-detection in a adjacent band of lower 
wavelength. In the case of the "J-band drop out" observed 
with HST, the object was observed in the 1.60 ji (H-band) 
but was not seen in the 1.15/.1 (Y -band) or 1.25 11 (J- band). 
The absence of emission in the latter bands is assumed to oc- 
cur due to Ly— a absorption by hydrogen clouds in between 
the source and us, allowing for an approximate estimate of 
the redshift of the object. More specifically we take as our 
dropout criterion 

AmAs ^ 1.2 (6) 

where Attiab is the difference in apparent magnitude 
between the two bands of observation, in this case the J and 
H bands. Observations at longer (near to mid IR) wave- 
lengths are required for photometric determination of ob- 
jects m ore distant than z = 1 0, ne cessitating JWST o bser- 
vations. iBouwens et al.l (120111 ) and lOesch et all (|201lh find 
a candidate z ~ 10 object in the co-added first and second 
year observation of the HUDF with the new WFC3/IR cam- 
era as a J-band dropout. This object is currently thought to 
be a galaxy, the most distant one observed to date, since 
the SED is a reasonable match to that of galaxies a t 2; > 9 
and it appears clearly extended (|Oesch et al.ll201ll ). How- 
ever the absence of spectra and the poor spatial resolution 
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Scenario Name 


Halo Mass Range 


Mns 


^form 


Zstart 


dn/dt 


dN 
dzdO'^ 




(Mo) 


(Mo) 






(Mpc-3yr-i ) 


arcmin~2 


A 


(1 - 2) X 10** 


10^ 


10 


13 


5 X 10"'-* 


235 


B 


(1 - 2) X 10** 


10^ 


12 


16 


7 X 10^1" 


16 


C 


(1 - 2) X 10** 


10^ 


15 


22 


1 X 10"^" 


0.77 


A 


(1 - 2) X lO'' 


10"^ 


10 


10.7 


5 X 10~** 


4435 


B 


(1 - 2) X lO'^ 


10^ 


12 


12.8 


6 X 10"** 


2965 


C 


(1 - 2) X 10^ 




15 


16 


2 X 10"** 


466 



Table 1. DM halo formation rates: dn/dt expressed in Mpc ^yr ^ , and J^^g-i as number formed per unit redshift and arcmin^ for 
cases considered in the text. Top three rows A-C are for a lO^'M© SMDS forming DM halos of mass (1-2) xIO^Mq ; bottom three rows 
are for the IO^Mq SMDS forming in DM halos of mass (1-2) xIO'^Mq. We have assumed that the DS started accreting baryons with a 
constant rate of lO~^M0/yr at Zstart and reached its final mass by Zform ■ 



of the image allow us to consider the possibility that this ob- 
ject could instead be a SMDS. Even though it may be hard 
to identify a DS uniquely with HST, the fact that at most 
one candidate has been found can be used to place bounds 
on the numbers of dark stars at redshifts up to z = 10. In 
this section we examine the observability of DS of various 
masses in existing HST imaging surveys, and in a later sec- 
tion examine the resulting bounds for future surveys with 
JWST. 



4.1 Comparison of DS stellar output with HST 
Sensitivity 



Figures l4l6| plot the predicted apparent magnitudes of Darli 
Stars of lO"* - lO^M© at various redshifts and compare these 
predictions to sensitivity of various HST surveys (plotted 
as thin horizontal lines) in two HST filters WFC3 F125 (J- 
band, colored blue) and F160 (H-band, colored red). In these 
figures, we have assumed that the SMDS formed at z=15 and 
survived to various redshifts as shown. In figures [4]and[5l the 
Dark Stars are considered to be formed via the extended adi- 
abatic contraction mechanism, without any captured DM; 
while in Figure [6] we consider the case with capture. 

The thick solid curves show the apparent magnitudes 
Mab for Dark Stars of various masses as a function of 
redshift in the J125 (F125W,blue) and Him (F160W,red). 
These solid curves are generated using simulated atmo- 
spheres spectra from TLUSTY (Fig [T|) and redshifting 
them, [Fu[\\ z)), imposing a cutoff at wavelengths lower 
than the Lyman-a if 2; > 6, assuming that photons 
at those wavelengths will be absorbed by the neutral 
hydrogen in the IGM. We use the H and J passbands 
throughput curves (T^-'(A)) for the HST WFC3, found at 
|http : //www. stsci . edu/-WFC3/UVIS/SystemTh roughput/, 
to compute the observed apparent magnitudes: 



-2.5 log 



j d\\T"'-'{\)F„{\;z) 
JdXXT"-J{\) 



-f-31.4 (7) 



The constant 31.4 is necessary to convert the fluxes to 
units of nJy. F{X; 2) is defined by: 



F,,{X;z) = 



{l + z)L,,{X') 
AnDliz) ' 



(8) 



where A is the redshifted wavelength, i.e. A = (1 -I- z)X' 
and L^, (A') is the emitted flux (we use TLUSTY to estimate 
it). The luminosity distance is labeled by Dl{z) and depends 



on the chosen cosmology. We define a J-band dropout to be 
any observation to the right of the green vertical line, cor- 
responding to a difference in apparent magnitudes of 1.2 or 
larger between the J and H filters as defined in Eq.(|B}, (the 
same criterion as used bv lOesch et al] (|201ll )). The location 
of the green line shows that J-band drop out technique will 
also identify the redshift of any SMDS found in this way to 
be at jz '--^ 10. 

In Figure |4] the sensi tivity limits f rom v arious deep 
field surveys compiled by lOesch et all (|201lf ): HUDF09, 
HUDF09-1, HUDF09-2, Early Release Science Data (ERS), 
CANDELS-Deep and CANDELS-Wide are indicated by dif- 
ferent line styles in the legends on the top right of each panel; 
these data are compared to the SMDS case of extended AC 
(no capture) . Also shown are the se nsitivity lim i ts for var- 
ious deep field surveys complied in lOesch et al.l (1201 ll ) . In 
Figure [3we focus on the most sensitive of these surveys, 
HUDFO9Q. 

Similarly, Figure |S] plots the apparent magnitudes as a 
function of redshift for lO'^M© (left) and IO^Mq (right) dark 
stars which grew via captured DM (rather than via extended 
AC). The SMDS formed with capture are harder to detect: 
since they are hotter their peak output is at lower wave- 
lengths (where Ly-a absorption is worse); in addition their 
radii are 5-10 times smaller, thus lowering their bolometric 
luminosities (|Freese et al.l [2010al ). In all cases the vertical 
dashed line is placed at the minimum redshift where the J 
band dropout criterion is satisfied. 

For SMDS with masses ^ 10^ Mq, the predicted fiuxes 
in both the F125W and F160W filters are too low to be 
seen in HST data; the only way around this would be if 
the objec t happened to b e gravi tationally lensed, as dis- 
cussed in IZackrisson et al.l (|2010aD . The 10^ Mq dark stars 
can be seen in the F125 (F160) passbands out to redshifts of 
9 (11.5) while the IO'^A/q dark star would be detectable out 
to redshifts of 10.5 ( 13). However, lO'^M© DSs would be too 
bright to be compatible with HST data: they would be sev- 
eral magnitudes brighter than the HST sensitivity, whereas 
the observed object is just bright enough to be seen. Thus 
the observed z — 10 candidate in HST cannot be a lO'^M© 
DS. In addition, if IO'^Mq SMDS formed at higher redshifts, 
we can place strong bounds on the numbers of them that can 
survive down to z=10, where they are not found. 



^ HUDF09 has a limiting 5a mab of 29.3 in the J-band for an 
exposure time of 94500 s and 29.4 in the H-band for an exposure 



time of 146711 s for the Hiqo band. 
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M|3s=l.e6Mo Extended AC Jjjj Dropout MDs=l-e7Mo Extended AC J, 35 Dropout 




Redshift Redshift 

Figure 4. SMDS in HST. Left (right) panels: Thick curves show apparent magnitudes in the H-band [F160W, sohd red] and J-band 
[F125W, blue curves]) for the 10"^ Mq (IO'^Mq) versus the redshift of observation for dark star formed via extended adiabatic contraction 
in a 10^M(ri (10 ^ M(7)) halo at redshift of 15. Thin horizontal lines indicate the 5(t detection limits of the various deep field surveys compiled 
bv lOesch et al ] l|201lh . with the areas of the surveys in arcmin indicated in the legends. The deepest survey to date is HUDF09 (lowest 
dotted horizontal lines). The vertical dashed line is placed at the minimum redshift where the J band dropout criterion is satisfied 
{z ~ 10). 



We also note that any SMDS that continued to exist to 
z=6 would have been seen as an 1775 dropout in HUDF which 
has a 29.9 niAB detection limit f or lOcr detection in the 1775 
passband ()Bouwens et al.|[2006l ). Since no candidates exist 
in the data, this makes it clear that SMDS did not survive 
to z=6. Thus we conclude that it is the l(f Mq SMDS that 
serve as the best possible explanation for the J-band dropout 
at z=10 seen by HST. 



4.2 Using HST observations to constrain the 
numbers of dark stars 

We will use HST data to constrain the fraction fsMDs(zstart) 
of early DM halos that can host SMDS. We focus on SMDS 
of masses Mds = 10^ — IO^Mq since lower mass DS are not 
observable in current HST data (unless they are significantly 
magnified by gravitational sensing or if they form clusters 

of dark stars ('Zackrisson et al."2010bjy). 

Following Zackrisson ct al. (2010bl. l20n"bl ). we compute 
the number Nobs of DS that could potentially be observed, 

Nobs ~ ^^^Q2 fSMDs{z — Zstart)0'^ fsurv fAt (9) 

and use the fact that at most one object has been observed 
with HST at z=10 to obtain bounds on fsMDs(zstart), the 
fraction of DM halos in a given mass range that can host a 
DS: 



Here dN/ dzdO^ is the number of DM halos forming per unit 
redshift per arcmin^ in which a given mass DS is hosted 



(computed from Figure We have multiplied by unit red- 
shift interval Az = 1, since we only consider SMDS formed 
within a redshift interval equal to one (see the discussion fol- 
lowing Eqn. (1)). Here 9^ is the total area surveyed in which 
the SMDS could have been detected, fsurv is the fraction 
of DS that survives from the redshift where the DS starts 
forming, Zstart , until it could be observed as a dropout (at 
2 ~ 10 with HST) and fAt is the fraction of the observa- 
tional window of time At during which the DS is still alive. 
Here, At is the cosmic time elapsed between the minimum 
and maximum redshift where the DS could be observed as a 
dropout. Please note that those redshifts are different from 
Zmin and Zmax defined under Equation [J] For the case of 
HST, we get At = 6.5 X lO'^yr (the cosmic time between 
the minimum redshift of 9.5 and maximum redshift of 10.5 
where the DS could be observed as a J-band dropout com- 
puted using Equation [5]) . 

We estimate the survey area 6^ in the following way: 
For each of the surveys in Figs. |4] and [G] we have indicated 
(in parentheses in the plots) the area (in arcmin^) observed 
by the survey. For DS of a given mass, we can add up the 
areas of all those surveys which are capable of observing DS 
as J-band dropouts to obtain a total effective area of ob- 
servability for that DS mass. In other words, we add the 
area of all surveys in which the fluxes in the //leo are still 
above the sensitivity limits while the fluxes in the J125 are 
a least 1.2 lower in apparent magnitude and below the de- 
tection limit of the J-band. From Figure |4] we estimate 
9^ — 4.7 X 3 arcmin^ as the effective area of the surveys 
in which a lO'^M© SMDS formed via extended AC could 
have been observed as a J band dropout with HST, since its 
only for the three deepest surveys, each with an area of 4.7 
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Figure 5. SMDS in HST: J125 (blue, F125W) and Hieo (red, F160W)) apparent magnitudes Mab for Dark Stars of mass ranging 
between ~ IO^Mq — IO^Mq as a function of redshift for WFC3 filters. Here the Dark Stars are considered to be formed via the extended 
adiabatic contraction mechanism, without any captured DM. The dashed horizontal lines represent the sensitivity limits for the deepest 
survey available, HUDF09. For the H band the Scr depth is 29.4 whereas for the J band it is 29.3. The exposure times are ~ 9.45 X 10^ 
s for the J125 field and ~ 1.47 X 10^ s for the Hieo field. The green vertical line corresponds to the lowest redshift where the dropout 
criterion is satisfied. Compared to Figure [4] now we explore a wider mass range for the SMDS. Note that SMDS of mass 10^ Mq or 
lower cannot be observed as J band dropouts with current HST data (another factor of 100 in observing time would be required) whereas 
heavier SMDS can be detected. 



arcmin'^ , that this SMDS would show up as a dropout. For 
the W^Mq SMDS formed via captured DM the detectabil- 
ity is much lower, implying that they could have been ob- 
served with HST WFC3 as a J band dropout only in the 
deepest survey, namely in HUDF09, which has an area of 
4.7 arcmin'^ . 

Although the z=10 J-band dropout seen by HST cannot 
be a 10^ M0 SMDS (as it would be too bright and would 



show up in both bands), still we can apply Eq. (|9| to place an 
upper bound the numbers of these objects. For the Mq 
stars formed via extended AC,, this area is increased to ~ 
160 arcmin^ , as all surveys compiled could pick this object 
up as a J band dropout. For the hotter DS fueled by captured 
DM, we can see from Figure [6] that the total area of the 
surveys in which lO^M© DS could have been detected is 
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Figure 6. SMDS in HST: Same as Figure |5]but for 10*' - lO^M© Dark Stars fueled by captured DM. The dashed horizontal lines 
represent the sensitivity limits for the deepest survey available, HUDF09. 



~ 160 arcmin^ (similar with the area for the extended AC 
DS of the same mass). 

We comment here on the three redshifts of formation we 
have chosen. For a conversion between Zform (redshift where 
the DS reaches its final mass) and Zstart (the redshift where 
the DS starts accreting baryons) see Table [T] 

• Case A: Zform = 10. Here, we assume that the DS be- 
come supermassive only at z=10 and not before. We can only 
constrain the product fsMDS x fsurv x /a*- The fraction of 
the observational window during which the DS is alive and 
can be observed, is /a* = min(r — rmin, At)/ At, where Tmin 
is the minimum DS lifetime that allows the DS to survive to 
z — 10.5 where it can be observed as a J-band dropout with 
HST. In the case of a IO'^Mq SMDS r„i„ ~ 1.15 x lO^rs 
(time elapsed between z — 13 and z — 10.5) whereas for the 
10*^ Mq SMDS 3.6 X 10 yrs (cosmic time elapsed 
between z = 10.7 and z = 10.5). We note that the limits we 
place on fsMDs(zatart) are only valid at Zstart ~ 13 (for the 
10'' M0 SMDS) and zstart ~ 11 (for the 10'' M© SMDS) as 
can be seen from Table [T] 

• Case B: zform = 12. Here we consider the DS to become 
supermassive at Zform ~ 12 and not at later redshifts. We will 
assume that the DS could survive until z ~ 10 {fsurv = 1) in 
order to constrain fsMDs(zstart) using null detection from 
HST Ji25 dropouts. From Table[T]we see that the Zstart value 
for the lO'^M© SMDS in this case is ~ 16 and for the 
lO'^M© SMDS it is ~ 13. In the case of a IO'Mq SMDS 
Tmin ~ 2.0 X 10**yrs whereas for the lOf' Mq SMDS Tmin ~ 
1.1 X 10**yrs. 

• Case C: Zfomi = 15. Here we assume the dark stars be- 
come supermassive by Zform ~ 15. The values for Zstart can 
be read off from Table [T] again. For the lO^M© SMDS 
Zstart ~ 22 and for the l<f Mq SMDS zstart ~ 16. This 
case is treated in a similar fashion as case B. For the 
10'' M© SMDS, Tmin = 2.9 X lO^yrs (the time elapsed be- 



tween redshifts 22 and 10.5) whereas for the lO'^M© SMDS 
"Train — 2.0 X 10 yrs (the time elapsed between redshifts 16 
and 10.5). 

From Eqn. (O, we obtain the following bounds for 
lO'M© SMDS formed via either extended AC or with cap- 
ture in each of the three cases (A-C): 

f-4.5 - \og{fsurv X /At), A 
-3.4-log(/,.u^„ X /At), B 
-2.1 " \og{faurv X /At), C 
(11) 

For lO'^M© SMDS formed via extended AC we get the 
following limits: 

f-4.8 - log(/s„ri, X /At), A 
-4.6-log(/,„,„ X /At), B 
-3.8 - \og{faurv X /At), C 

(12) 

The values of Zstart that correspond to these values of 
Zform can be found in the last three rows of Table [1] The 
reason that the bounds on the numbers of lO'^M© SMDS are 
tighter than those on the lO'M© SMDS is the following. In 
order to reach a larger mass, the DS had to start forming 
at an earlier redshift and in larger halos; but the numbers 
of larger halos that can host DS is smaller at higher red- 
shifts. Similarly, the bounds in Case A are ~ 10 (~ 300) 
times stronger than the bounds in Case C for the 10® M© 
(IO'M© ) SMDS. Again the reason for the very large dis- 
crepancy ~ 300 is the fast decrease of the formation rate of 
1 - 2xlO*M© DM halos at redshifts higher than 2: ~ 15 as 
can be seen from Figure [1 For SMDS lighter than lO'^M© 
HST cannot be used to place constraints, as those objects 
are not detectable with HST as J-Band dropouts. A sum- 
mary of our bounds can be found in Figure [7] where we plot 
the exclusion limits for fsMDS • 
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Figure 7. Upper bounds from HST on the fraction fsMDS of early halo hosting dark stars of masses as labeled above the plots; values 
above the lines are excluded. Different lines correspond to different values of the redshift (zfoj-m = 10, f 2, and 15) at which the DS attains 
this mass; see Table [T]for the connection between (a) the redshift Zgtart at which the DS came into existence and started to grow and 
(b) the redshift zjarm a-t which it reached the supermassive size as labeled. This plot assumes fsurv =1 (i-e. the DS survives long enough 
to reach the redshift window of observability as a J125 dropout with IfST). However, the DS need not survive throughout the entire 
window; in fact the horizontal axis in both plots is logj^Q /a*, for which a value of corresponds to the DS lifetime being sufficiently 
large that it survives throughout the redshift window of observability. Solid lines correspond to DS formation via extended AC (without 
capture) while dashed lines correspond to DS formation via capture. Since DS less massive than 10^ Mq are too faint to be detected by 
HST, these data do not bound fsMDS for lower mass dark stars. 



4.3 Other Bounds on Numbers of SMDS 

Further bounds on the numbers of DS and the halos they 
form in should result from a variety of considerations. 
One would be th e contribution to reionization. Work of 
IVenkatesanI (|200Cll ) studied stellar reionization with the stan- 
dard fusion powered first stars (Population III), without any 
dark stars. From comparison with the optical depth to last 
scattering from early WMAP data, she bounded the fraction 
of baryons in halos that can cool and form stars (assuming a 
Scalo initial mass function) to be in the range /* ~ 0.01—0.1. 
However, it is not clear how these numbers would change in 
the presence of DS and with the updated value f or the op- 
tical depth from WPMAP7 (jKomatsu et al.ll2011^ ■ 

More recently, the effects of DS (and the resul- 
tant main sequence stars) o n re ionization was stu died by 
ISchleicher et~ai] (|2008l . |2009| ) and lScott etHI (|201ll ). While 
DS are fully DM powered, they remain so puffy and cool 
that no ionizing photons are produced, and there is no con- 
tribution to reionization. However, once the DM fuel begins 
to run out, they contract and heat up as they approach the 
zero age main sequence (ZAMS) with the onset of fusion, 
at which point they do produce ionizing photons. For the 
case of extended A C, and for DS less massive than lOOOAf©, 
IScott et al.l (|201ll ) concluded that the reionization history of 
the Universe is unaffected by the DS, compared to the case 
of more standard Pop HI stars: the DS period of no ion- 
izing photon production is compensated by a short period 
of high ionizing photon production during approach to the 
ZAMS. However, we are not sure what the effect on reion- 
ization would be in the case of the more massive SMDS. On 
the one hand the more massive stars are hotter and brighter 
and would emit substantial amounts of ionizing photons; on 



the other hand the more massive the star, the shorter the 
lifetime. 

Fo r the case of DS with high capture rates, previous 
studies lScott et al] (|201ll ) find that reionization is somewhat 
delayed, decreasing the integrated optical depth to the sur- 
face of last scattering of the CMB. However, variation of 
astrophysical parameters for the case of standard reioniza- 
tion with standard Pop HI stars can produce exactly the 
same effect, so that disentangling the se effects will prove 
difficult. Nonetheless IScott et al.l l|201lh do argue that they 
can rule out the section of parameter space where dark stars 
~ IOOOM0 with high scattering-induced capture rates tie 
up more than 90% of all the first star-forming baryons and 
live for more than 250 Myr. Again, their work should be 
extended to the heavier SMDS we study in this paper. 

A complicating factor (for both the cases of extended 
AC and capture) is that the SMDS do eventually col- 
lapse to BH, and it's not clear how rapidly that happens. 
If the collapse to BH is rapid, this may cut short the 
ZAMS phase and reduce the role SMDS play in reioniza- 
tion. Second, the SMDS are likely to have stellar pulsations 
(jMontgomerv et al.ll201ll . in progress); as a consequence it is 
possible they will lose some mass before reaching the ZAMS. 
Third, even after joining the ZAMS, en route t o BH collapse, 
the SM DS may blow off some of their mass (|Umeda et alj 
(j2009l) suggest 1/2 of their mass). 

Alex Heger (private communication) has the following 
new results for early stars (only made of hydrogen and he- 
lium) that are nonrotating; If they are heavier than 153,000 
Mq, no hydrostatic equilibrium solution exists, i.e. no pri- 
mordial hydrogen burning star exists. Thus once a fusion 
powered star accretes enough mass to get heavier than this, 
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then it collapses straight to a BH. For any of our dark stars 
that are heavier than this, once they run out of DM, they 
collapse directly to BH without contributing at all to reion- 
ization. Rotation might change these results. 

Further, there are implications of DS regarding the frac- 
tion of baryons that end up in DSs. Our work assumes that 
the DS can grow in a DM halo of a given mass until almost 
al the baryons in the halo (assumed to be the baryonic mass 
fraction in the Universe) are accreted onto the DS. If the to- 
tal fraction of halos in which such DSs form is too large, this 
implies that most of the baryons in the Universe are trapped 
inside DS and and it is not clear how they would contribute 
any further to galaxy formation. As mentioned above, en 
route to BEI collapse, the SMDS may blow off some of their 
mass, reinjecting baryons into the surrounding halos and 
alleviating this problem somewhat. 

Furth er bounds on the nurn bers of dark stars have been 
studied in ISandick et all l|201ll ) The remnant black holes 
from the DS should still exist today, including inside the 
Milky Way. They still have enhanced amounts of DM around 
them, known as DM spikes. The DM inside the spikes an- 
nihilates to a variety of final products, with 7-rays that 
would be detected by the Fermi Gamm a Ray Space Tele- 
scope (FGST). In ISandick et al.l |2oi3), it was noted that 
most of the 368 point sources observed by FGST might 
in fact be due to DM annihilation in the spikes. In addi- 
tion, FGST data were used to place bounds on the frac- 
tion of early haloes hosting DS to avoid overproduction of 
gamma-rays from annihilation in the remnant DM spikes. 
The bounds range from fos < 10^'^ — 1, depending on the 
WIMP mass and annihilation channel. 

All of these considerations are beyond the scope of this 
paper. For now we take these arguments to imply that not 
every early halo can contain a DS. 



JWST NIRCam filters. Figures [8]and[9]plot the stellar spec- 
tra of SMDSs of various masses and formation redshifts as 
a function of wavelength (for light emitted at z = 15, 10 
and 5) and compare to the sensitivity of JWST filters for 
lO^'s and lO^s exposure times. In Figures 1101 and 1111 we 
instead plot the apparent magnitudes as a function of red- 
shift of emitted light for various SMDS through the NIR 
camera wide passband filters, with each panel in the figure 
focusing on a particular JWST broadband filter; in these 
two figures the SMDS are formed via extended AC and cap- 
ture respectively. Ly— q absorption cuts off the photons with 
wavelengths lower than 1216A(in the rest frame); we treat 
the absorption as being complete. Thus the SMDSs drops 
below the JWST sensitivity limit at 2 ~ 6 for the F070W 
filter and at z ~ 10 for the F115W case. 

Since the most massive dark stars are the brightest, they 
are the easiest to detect. From figures 1 101 and [TT1 one can see 
that 10^1^0 dark stars, both with and without capture, are 
individually observable in 10^ seconds of NIRCam data even 
at redshifts as high as 15 in filters with a passband centered 
at 2 ^m and higher (F200W-F444W filters). For the case of 
a Mq SMDS, a longer exposure time of lO^s allows the 
dark star, both with and without capture, to be individually 
observable in all filters from F200W to F444W even sk, z 
15. For IQ^Mq SMDS, those formed via extended AC are 
visible in these filters out to 2: ~ 15 with 10^ sec exposure 
time while those formed with capture are too dim. Lighter 
ones < 10^ M0 would not be detectable as J-band dropouts 
but, if they survived to lower redshifts (e.g. z=7) would likely 
already have been seen with HST or other telescopes. Since 
the sensitivity of the higher wavelength MIRI filters above 
5/i is worse, only the 10^ Af© DS are bright enough to be 
observable in MIRI filters (see the discussion of Figure [T^ in 
Section |6]). 



5 OBSERVING SUPERMASSIVE DARK 
STARS WITH JWST 

Dark stars can be detected by upcoming James Webb Space 
Telescope (JWST). Table [2] gives a summary of the sensitiv- 
ity of the NIRCam and MIRI cameras on JWST in various 
wavelength bands 0. One can see that the NIRCam is much 
more sensitive than the MIRI filters, so that light emitted 
at wavelengths larger than 5 microns is harder to observe. 
In this section we estimate the number of SMDS that would 
show up in a typical survey with JWST NIRCam or MIRI 
cameras, based on the bounds we have just derived in the 
previous section. 

Figures [Sl-fTTI illustrate the detectability of SMDS with 



^ Specifically we show the lOcr required niAB sensitivities 

for the NIRCam and MIRI wide filters after 10^ seconds 

exposure derived based on the limiting fluxes published at 

http : //www. stsci . edu/jwst/ instrijjnents/iniri/ sensitivity/ 

and http: / /www. stsci . edu/jwst/ ins truments/nir cam/ sensitivit 

One can scale the limits to different exposure times, as the lim- 
— 1/2 

iting flux (X tgxpo3ure and converting to rriAB magnitudes is just 
a matter of applying Equation [7] For instance an increase by a 
factor of 100 in exposure times would convert in a gain in the 
sensitivity limits by 2.5 AB magnitudes. 



5.1 Detection at 2 ~ 10 as a J115 band dropout 
with JWST 

As expected, the dark stars that could have been detected 
as J-band dropouts with HST are also detectable using the 
same technique with JWST. We again will adopt the same 
dropout criterion, i.e. AniAB > 1-2 in the J and H broad- 
band filters. Figure [T^] shows the sensitivity of JWST in 
a 10'* exposure in the 1.15^ (J-band) and l.bO/i (H-band) 
filters for NIRCam. The apparent magnitudes for 10^ and 
IO^Mq SMDS with and without capture are also shown 
for comparison. Here, the SMDS form at z=15 and are as- 
sumed to survive to various redshifts as shown. Compar- 
ing Figures |4] and [12] (see also Table [2]) , one can see that 
JWST is about half a magnitude more sensitive than HST 
to finding SMDS as J-band dropouts (for 10"'s exposure time 
with numbers provided in the literature as 5(7 detection in 
HUDF09 and lOa detection with JWST). 

The three cases of IO'^Mq with or without capture as 
well as 10^ Mq without capture could be detectable in a 
JWST survey as J band dropouts in the redshift range 
A>^e* ^'^ even with the lower 10* second exposure times. 
The 10^ M0 Dark Star formed via captured DM (lower left 
plot) in Figure [12] will require a longer exposure time of 
lO'^s (which would correspond to the same exposure time as 
the 2004 HUDF survey). 

In order to predict how many SMDS would be visible 
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Filter 


^center { ) 


lo§lQ ^center 


AA ( nm) 




NIRCam 


F070W 


0.7 


-0.15 


0.175 


28.1 




F090W 


0.9 


-0.05 


0.225 


28.51 




F115W 


1.15 


0.06 


0.2875 


28.72 




F150W 


1.5 


0.17 


0.375 


28.77 




F200W 


2.0 


0.30 


0.5 


28.75 




F277W 


2.77 


0.44 


0.6925 


28.67 




F356W 


3.56 


0.55 


0.89 


28.55 




F444W 


4.44 


0.65 


1.11 


27.92 


MIRI 


F560W 


5.6 


0.75 


1.2 


25.65 




F770W 


7.7 


0.89 


2.2 


25.28 




FIOOOW 


10.0 


1.00 


2.0 


24.29 




F1130W 


11.3 


1.05 


0.7 


23.32 




F1280W 


12.8 


1.1 


2.4 


23.53 




F1500W 


15.0 


1.18 


3.0 


23.26 




F1800W 


18.0 


1.25 


3.0 


22.31 




F2100W 


21.0 


1.32 


5.0 


21.56 




F2550W 


25.5 


1.4 


4.0 


20.28 



Table 2. lOcr sensitivity limits for the JWST wide passband filters (fourth column). The mjiB limits are derived assuming 10* sec- 
onds exposure and are based on the limiting fliaxes published at ^http: //www. stsci . edu/jwst/instrinnents/miri/sensitivity/ and 
|http : //www. stscl . edu/jwst/ lnstruments/nircam/sensltivity/table*! We identify each filter by its name, in the first column. Values 
in the second column correspond to the center wavelength of each filter, whereas in the fourth column the values for the passband width 
are given. 



in a JWST deep field survey we liave to assume sometliing 
about tlie total field of view (FOV) of all future JWST 
surveys in which the stars would be observable. The FOV 
of the NIRCam instrument is 2.2' x 4.4' = 9.68 arcmin^ (see 
]http : //www. stsci . edu/jwst /overview/design/) This 
value is likely to be an underestimate. Since HST had 
multiple surveys with a total of 160 arcmin^ , we will also 
consider the case of multiple surveys with JWST with a 
total added area of ~ 150 arcmin'^ . Given the bounds on 
the numbers of DS from HST from the previous section, 
we find that the number of expected SMDS with JWST 
as J-band dropouts is N < 1 and therefore conclude that 
SMDS are hard to detect with JWST as J-band dropouts. 
This is expected since HST was already sensitive enough to 
observe them as J-band dropouts, assuming enough would 
have survived from their formation redshift until z ~ 10. 
The only improvement could be made by a larger survey 
area compared to the one with HST. For the IO^Mq 
SMDSs formed with capture, which were detectable only 
in the 4.7 arcmin^ of HUDF09, JWST should be able to 
provide a larger survey area so that these objects become 
more detectable. 



5.2 Detection at z ~ 12 as a Hiso band dropout 
with JWST 

Whereas JWST is not particularly better than HST at find- 
ing J-band dropouts, it will be significantly better at find- 
ing SMDS as H-band and K-band dropouts at higher wave- 
lengths. In this section we focus on H-band dropouts, where 
the object can be seen in the F200 NIRCam fiher of JWST 
but not in the F150 NIRCam fiher. As before in Eqn. 
we require the difference between the broadband fluxes in 
the Hi5o and K200 filters to be greater than 1.2 AB magni- 
tudes. We see that the SMDS stellar light seen with JWST's 
if 200 filter is essentially unaffected by Ly— a absorption un- 
til z ^ 15, whereas the IGM absorption will cut off most 



of the fiux in the H150 at 2; > 11.5. (see Figure [TS}. Hence 
SDMS can appear as H-band dropouts. 

We will consider the case of SMDS forming at zform =12, 
the same as the time of observation. Figure [13] shows that 
the three cases of lO^M© SDMSs with and without cap- 
ture as well as 10^ M© SMDSs without capture are all de- 
tectable in a JWST survey as H150 dropouts in the redshift 
range 11.5 — 12.5. DS formed at higher redshifts could be 
seen all the way out to z ~ 14, 15, but likelihood analyses 
on any objects found as H-band dropouts with photometry 
with JWST will probably estimate the redshift at z '-^ 12. 
The 10^ M0 Dark Star formed via captured DM (lower left 
plot) is too faint to appear as a dropout. The number of 
Hi5o dropout events is given by Eq. ((2} with Az = 1 and 
fsurv = 1 since the objects are observed at the same time 
they form) and the appropriate survey area 6^ for JWST 
must be applied. 

Is it reasonable to apply the bounds from HST on the 
numbers of SMDS at z=10 to those at z = 12? We will con- 
sider three different possibilities, and summarize all results 
for the predicted number of H band dropouts with JWST in 
Table |3l If we assume that all the SMDS at z = 12 have the 
same properties as those at z — 10, and that they survive 
throughout the redshift window observable by HST, then 
the HST bounds are so stringent that JWST will not be 
able to see many of them. This is the case we label "Maxi- 
mal Bounds" . In particular, IO^Mq SMDSs would have been 
so bright as to be easily seen in HST, and the resultant strin- 
gent bounds imply that only Nobs ~ 1 DS would be found 
even with multiple surveys with 150 arcmin^ FOV. For 
IO'^Mq SMDSs the bounds from HST are shghtly weaker 
because the objects are not as bright, so that 10 (32) of 
these might be found per 150 arcmin^ field for DS that grew 
via extended AC (with capture). Since the ones with capture 
are fainter and harder to see, (counterintuitively) the weaker 
HST bounds imply that more of them might be found with 
JWST. 
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Figure 8. Spectra for supermassive dark stars formed at Zform = 15 compared with sensitivity of JWST filters. Listed above each panel 
are the mass of the DS in solar masses, the formation mechanism (extended AC or "with capture") and the surface temperature T^ff. 
The fluxes are shown at z = lb (dashed line), 10 (solid line) and 5 (dotted line) and compared to the detection limits of NirCam wide 
passband filters. The colored horizontal lines represent the sensitivity limits for the filters as labeled in the legend for exposure times 
10*s (upper lines) and lO^s (lower lines). IGM absorption will decrease the observed fluxes for wavelengths shortward of the vertical red 
lines, which indicate the Ly-a line (1216A) redshifted from the rest-frame of the star. 



However, it is very likely that there are more SDMS at 
z=12 (the JWST window) than at z=10 (the HST window). 
For one thing, the host halo formation in this mass range 
peaks at 2 ~ 12 (see Fig 3). Moreover at lower redshifts 
{z ~ 10) the DM halos that could host those SMDSs are 
much more likely to merge to form even larger halos. In ad- 
dition, after the first SDMS die (before z=10), they turn into 
fusion powered stars that produce ionizing photons, which 
disrupt the formation of DS at lower redshifts. Indeed the 
strong halo clustering at high redshift would cause the pos- 
sible formation sites to be preferentially close to or inside 
the HII regions during reionization, potentially leading to 
strong sup pression of star formation; due to this mechanism 
llliev et al.l (fibo?, ) found a suppression of 10* Mq -10^ M© 
halos by an order of magnitude due to Jeans mass filtering 
in the ionized and heated H II regions. 



We will thus recalculate the number of DSs detectable 
with JWST using weakened bounds from HST. We will take 
fsMDS fAt fsurv = 1.5 X 10~^ as our "Intermediate Bounds" 
case. This case could imply that not all minihaloes host DS, 
or that not all DS survive throughout the z~10 HST observ- 
ability window. In this Intermediate Bounds case, hundreds 
or thousands of SMDS are potentially observable. 

For comparison, in the Table we list as a third case 
the full number of DM haloes that could in principle host 
DS. If all of these contained DS one would expect up to 
~ 450, 000 DS with JWST. However, as discussed in Section 
14.31 this would be extremely unreasonable as there would 
be no baryons outside of DS left for galaxy formation. Our 
results for the detectability of SMDS as H-band dropouts 
with JWST are summarized in Table [S] 
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Figure 9. Similar to Figure |8] now for a 10^ Mq DS formed eitiier at Zf^ 
a 10^ M0 DM iialo (rigiit panel). 



=20 in a IO'^Mq DM halo (left panel) or at 
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Table 3. Upper limits on the number of SMDS detections as -f/150 dropouts with JWST. In first three rows (labeled "Maximal Bounds") 
we assume that all the DS live to below z=10 where they would be observable by HST, and we apply the bounds on the numbers of DS 
fsMDS from HST data in Section 14.21 The middle three rows (labeled "Intermediate") relax those bounds by assuming that only ~ 10~^ 
of the possible DS forming in z=12 haloes make it through the HST observability window. For comparison we also tabulate in the last 
three rows the total number of potential DM host halos in each case. We also split the number of observations in two categories, N^^^ 
and A^^"'**. The first assumes a sliver with the area equal to the FOV of the instrument (9.68 arcmin^ ), whereas in the second we 
assume multiple surveys with a total area of 150 arcmin'^ . Note that for the case of the IO^Mq SMDS the predictions are insensitive 
to the formation mechanism. 



5.3 Detection at z ~ 15 as a 7^200 dropout with 
JWST 

DS at z > 14 can be detected as ^"200 band dropouts using 
the F200 and F277 NirCam filters in JWST, as shown in Fig- 
ure [Hfor 10^ and IO'^Mq SMDS formed via extended AC 
(no capture) at zform =20. To qualify as a ^^200 dropout the 
difference in magnitudes between the F277W and F200W 
filters must be greater than 1.2. As for the case of H-band 
dropouts above, we use HST data to bound the number 
of possible K-band dropouts, under three different assump- 
tions: (i) Maximal Bounds where every DS survives through 
the HST observability window at z~10; (ii) Intermediate 
Bounds with ~ 10~^ of the possible DS surviving that long, 
and (iii) for comparison simply counting every possible host- 
ing halo. Our results for predicted numbers of SMDS observ- 
able as K-band dropouts with JWST are summarized in in 
Table |4] 

The 10® Mq DS could be observed in the redshift range 



z ~ 15 ~ 17 as a 7^200 dropout for 10* seconds exposure. For 
the case of Maximal Bounds from HST, we predict at most 
A^Sr"' ~ 1- For the Intermediate Bounds case, the possible 
number of detections is increased to roughly 5 for the case of 
a 9.68 arcmin^ FOV or to 75 for the 150 arcmin'^ case. The 
(unreasonable) case where every possible halo hosts a DS 
shows the maximal number of 10® Mq SMDS observable as 
K-band dropouts to be ~ 70, 000. In the case of the IQ'Mq 
star, it would appear as a ^^200 dropout in the 16 — 20 red- 
shift range. However due to the sharp drop in the formation 
rate of DM halos in the 1 — 2x lO^M© at such high redshift 
the number of dropout events we predict in this case is at 
most ^ 1 (other than for the unreasonable case where ev- 
ery single possible halo hosts a SMDS). The results for the 
detectability of SMDS as K-band dropouts with JWST are 
summarized in Table |4] 
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Figure 10. Apparent magnitudes as a function of redshift for various SMDS through the NIR camera wide passband filters on JWST 
for the case of formation via extended AC. The number after the letter F and before the letter W in the name of each filter corresponds 
to the wavelength in the center of the passband in 0.01 ytm units. The two horizontal lines correspond to sensitivity limits for each filter 
for lO^s exposure time (the dotted line) and 10^ s exposure time (the dash-dotted line). The zform labeled in the legend is the formation 
redshift when the SMDS reached its corresponding mass. The curves corresponding to Zfo^m = 15 do not extend all the way to z= 20 
because at that high redshift the star has not formed yet. The sharp decrease of the fluxes at various redshifts in the first three panels 
is due to the Gunn-Peterson trough entering the filters. The higher wavelength filters F277W-F444W would not be affected by the IGM 
absorption until z > 20. 
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Figure 11. Same as Figure [Tol for SMDS formed "with capture" in various JWST bands as labeled. 
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Figure 12. SMDS with JWST as J115 band dropouts; Apparent magnitudes for various SMDS tlirougli ttie F115W and F150W filters 
for NirCam. Top panel: 10^ Mq and W Mq Dark St ars formed witliout DM capture. Lower panel: 10^ Mq and IO'^Mq Dark Stars formed 
"with capture". The dotted horizontal lines are obtained from the lOcr required sensitivities for 10'* seconds exposure data published at 
|http: //www. stsci ■ edu/jwst/ instruments/nircam/ note that the detection limits for the J115 and the -ffiso filters differ by only ~ 0.05 
rriAB apparent magnitude and arc thus essentially indistinguishable. The dashed horizontal lines are obtained assuming 10^ seconds 
exposure time. 
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Table 4. The number of SMDS detections as K200 dropouts with JWST. Cases are the same as above in Table [3] 



6 SMDS VS POP III GALAXIES WITH JWST 

A key question in the discovery of dark stars with JWST 
will be the ability to differentiate these objects from other 
sources at high redshifts. Assuming that a population of 
potential 2 > 10 candidates is identified by the drop out 



techniques described in previous sections, the most signifi- 
cant contaminant population at these redshifts is likely to 
be galaxies dominated by Pop III stars. In this section we 
focus on ways to different iate between SMDSs and galaxies 
containing Pop III stars. IZackrisson et al] (l2010al ) showed 
that DSs in the mass range < IQ^ Mq could be easily distin- 
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Figure 13. SMDS with JWST as Hi^o band dropouts: Apparent magnitudes for SMDS through the F150W and F200W NirCam filters. 
Those could be used to establish dropout detection criteria in the 12 — 14 redshift range. Top panel: cases of interest (10"^ and IO'^Mq) 
Dark Stars formed without considering DM capture. Lower panel: 10® Mq and 10^ Mq Dark Stars formed including DM capture. The 
vertical green dashed line indicates the minimum redshift at which the DS will appear as a dropout. 



guished from galaxies in the redshift range 2; = — 15 (in- 
cluding galaxies containing Pop III stars), supernovae, AGN, 
Milky Way halo stars as well as Milky Way brown dwarfs 
by their extremely red colors in color-color plots. The DS 
considered there have all Te// < 9000A', which leads to an 
decrease of the ratio (B) of the fluxes to the left and right 
of the Balmer j ump located at 0.365 [im with temperature 
(see Sec. 8.3 of iRuttenllioO^ V The significant Balmer jump 
in the case of DS with mass <10^M© will lead to very red 
JTiaes — 771444 c olors at z=10, offe r ing a d istinct signature, as 
pointed out in IZackrisson et "all (|2010al ). Here we study in- 
stead the much heavier SMDSs with Mds > 10^ Mq. These 
heavier stars are intrinsically much brighter and thus eas- 
ier to find as dropouts. Ifowever, they are also hotter than 
lOOOOif, leading to values of the ratio B to increase with 
temperature, as explained in Sec. 8.3 of iRuttenI (|2003l ). For 
the SMDS we consider here the Balmer jump is insignifi- 
cant, therefore its much more difficult distinguishing them 
from potential interlopers based on the technique proposed 



in'Zackris son et al.l l|2010ah for the smaller ~10^M© DS. In 
this section we begin a discussion of this issue, restricting our 
studies to what can be learned from JWST directly. Future 
studies will be required, in which we investigate also the pos- 
sible role of spectroscopy with TMT and GMT, and other 
upcoming observatories in differentiating Pop III galaxies 
from dark stars. 

The earliest Population III stars (in the absence of 
dark matter heating) are expected to have masses in the 
range 10-100 Mq - too faint to be seen as individual ob- 
jects with JWST llqh|[l999l : IOh" et al. I I2OOII : iGardnCT et al.l 
I2OO6I : iRvdberg et al.l 2010l ). However a galaxy containing 
10° - 10'^ Mp) of Pop III stars might indeed be detectable. 
IZackrisson et al] l|2011ah presented a comprehensive study 
of the integrated spectra signatures of Pop III stars in the 
wide filters of JWST. Their main findings are that Pop III 
galaxies could be detectable to redshifts as high as 20 if 
the stellar population mass is ~ 10^ Mq (or in the case of 
10° M© stellar population mass up to redshifts of 10). A 
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Figure 14. SMDS with JWST as ^200 band dropouts: Apparent magnitudes for SMDS formed without DM capture through the F200W 
and F277W NirCam filters. Left panel: for the IO^A/q dark star. Right panel: for the lO' Mq dark star. The vertical green dashed line 
indicates the minimum redshift at which the DS will appear as a dropout. 



similar study by iPawlik et al.l (|201ll ') , who examined neb- 
ular emission lines from early galaxies, came to the same 
conclusion : thousands of the se m ay be found with JWST . 
Moreover, llnoue et all (|201l| ) and IZackrisson et al.l (|2011al l 
have pr oposed selec t ion c riteria using two of the filters of 
JWST: llnoue et al.l (1201 ih argued for using two NIRCAM 
filters and Zackrisson et al.l ()2011al l argued for adding imag- 
ing in two MIRI filters to more cleanly differentiate between 
Pop III galaxies and Pop II or Pop I galaxies at z~ 7 — 8. 
Indeed these authors found that galaxies containing Pop III 
stars at high redshift are typically brighter in most JWST 
filters than later generations of stars; thus galaxies with Pop 
III stars would be the most likely source of confusion in iden- 
tifying dark stars. 

Using the Yggdrasijf] model grids (jZackrisson et al.l 

l2011al . see http : / /ttt . astro . su. se/~ez/[) we compare sig- 
natures in the NIRCam passbands of Pop III galaxies at 
z~ 10 — 15 with those of SMDS. All the nomenclature use d 
here for Pop III galaxies follows IZackrisson et al] (|2011al V 
we consider three different Initial Mass Functions (IMF) for 
Pop III galaxies: 

• Pop III.l: A zero-metallicity population with an ex- 
tremely top hea vy IMF and a Single Stellar Popula- 
tion (SSP) from ISchaereJ (|2002D with a power-law IMF 
{dN/dM oc M~°'). The population has stellar masses in 
the range 50 — 5OOM0 and a Salpeter slope a — 2.35 for the 
entire mas range. 

• Pop III. 2: A zero-metallicity po pulation wit h a m od- 
erately top-heavy IMF. A SSP from iRaiter et all (|2010l ') is 
used. This model has a log-normal IMF with characteristic 
mass Mc — lOM© and dispersion a = IMq. The wings of 
the mass function extend from 1 to SOOM©. 

• Pop III, Kro upa IMF: In view of recent simulations (e.g. 
iGreif et aklboid ) the mass of Pop III stars might be lower 

* We highly recommend watching the movie Thor to understand 
this name. 



than pr e vious ly predicted. Therefore in this case a normal 
iKroupal |20o3) IMF, usually describing Pop II/I galaxies, is 
used. The stellar masses range in the 0.1 — lOOMp and the 
SSP is a rescaled version of the one used in ISchaereil (|2002l ) 

Following IZackrisson et al] (|2011al ') we further subdi- 
vide the models into two types, based on the amount of 
nebular emission. The first galaxies are expected to have 
significant ionized gas surrounding them. Depending on how 
compact the HII region is, the escape fraction for ionizing 
radiation from the galaxy into the IGM can vary anywhere 
from 0-1. Hence we consider the two extreme possibilities: 

• Type A galaxies: If the gas covering fraction fcov ~ 1, 
then there is maximal nebular contribution to the SED and 
no escape of Lyman continuum photons. 

• Type C galaxies: If fcov ~ 0, there is no nebular con- 
tribution to the SEDs and instead stellar light dominates 
the SED. We will not consider here the intermediate case of 
Type B galax;ies. 

IZackrisson et al] (l2011al ) argue that the nebular emis- 
sion typically dominates the spectrum of young Pop III 
galaxies at z~ 10 — 15; e.g. at z=10 nebular emission dom- 
inates for galaxies younger than 10 Myr. All young or star 
forming galaxies are expected to have significant contribu- 
tion to the SEDs from nebular emission, and this effect is 
increased with lower metallicity or a more pronounced top 
heavy IMF. Hence we will predominantly focus on Case A 
of maximal nebular emission from the early galaxies. 

In Figure [15] we plot the SEDs (in apparent magni- 
tudes) of SMDS and Pop HI galaxies at z=12 as a function 
wavelength. Our interest is in their detectability with the 
NIRCam and MIRI cameras on JWST. The vertical dotted 
line demarcates the wavelength ranges covered by the two in- 
struments, and the dark blue horizontal segments represent 
band widths and the sensitivity limits of individual filters 
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Figure 15. JWST detection limits and apparent magnitudes of SMDS and Pop III galaxies. Each filter's sensitivity limit for JWST 
NIRCam and MIRI wide filters for 100 hours of exposure time is plotted as a segment at the corresponding wavelength. Left panel: 
Apparent magnitudes for SMDS in the IO^Mq (solid lines)-lO'^Af0 (dash-dotted lines) mass range for both extended AC (red lines) 
and Capture (green lines) formation mechanisms. Right panel: Apparent magnitudes of SMDS of 10^ Mq (solid lines) and Pop III 
instantaneous burst galaxies of 1 Myr age (dashed lines) with the same stellar population mass as the SMDS. 



assuming a 100 hour exposure 0. In the left panel we plot 
the apparent magnitudes for 10^ and IO^Mq DS formed via 
both Extended AC and Capture mechanisms. We have pre- 
viously discussed (see Figures 

[lO]and[Tl| that both 10^ and 
W'^Mq DS are bright enough to be observed by the NIR- 
Cam filters. On the other hand, in the less sensitive MIRI 
filters, IO'^Mq DS can be seen in the lowest two wavelength 
filters (F560W and F770W) but lO'^M© DS are too faint to 
be observed. 

In the right panel we compare the observed SEDs of 
lO** Mq SMDS (solid curves) with Pop III galaxies (dashed 
curves). For the galaxies we assumed an instantaneous star- 
burst (at t = 0) and used the results from the Yddgrasil 
code at 1 Myr after the burst. The light from the galaxies 
is assumed to be dominated by nebular emission (Type A) 
l|Zackrisson et all boila) for galaxies younger than 10 Myr. 
We have taken the stellar mass of the galaxy to be the same 
value 10^ M0 as the DS mass. SMDS are brighter than the 
galaxies in all filters in which the objects are potentially vis- 
ible. The sharp cut off in flux at logjQ Xobs ~ 0.02 is due to 
Ly-a absorption. For a stellar population mass of W'^Mq , 
Pop III.l galaxies are detectable as a i?i5o dropouts in a 
deep field survey with an exposure of 100 hours; Pop III. 2 
are still just above the sensitivity limits; but Pop III galaxies 
with a Kroupa IMF are not detectable as -ffiso dropouts. 

Let us imagine that an object has been detected as 
a photometric dropout at some redshift, say an H-band 
dropout at z — 12. Our goal is to determine the nature 
of this object, i.e., to differentiate SMDS from first galaxies 
with JWST. One approach would be to exploit the emission 
lines i n galaxies that are not shared by the DS. IPawlik et al] 
l|201ll ) have shown that there would be several major sig- 



5 For NIRCam we did not plot the F090W filter, since the 
throughput profile is not yet available. 



natures in the spectra for Pop III galaxies with significant 
nebular emission (our Case A), including the Hell line at 
0.1640 [im and Ha emission. They found that JWST spec- 
trometers (NIRSpec and MIRI) are indeed sensitive enough 
to detect these emission lines, thereby potentially finding up 
to tens of thousands of star-bursting galaxies with redshifts z 
>10 in its field of view of ~ lOarcmin^. They also found that 
the Hel640 recombination line is only detectable in signifi- 
cant numbers for the case of zero-metallicity starbursts with 
top- heavy IMF. They noted that their estimates are consis- 
tent with previous est i mates of JWST starburst counts (e.g 
iHaiman fc Loeb|[l99^ : [Ohlll999l ). A third possibility would 
be to detect the continuum limit of the Balmer series at 
0.3646 \im in the rest frame. 

In short, if foUowup spectroscopy is done on an object 
found as as dropout with JWST, the detection of a Hell 
1640 emission line or an Ha emission line would most likely 
indicate that the object is a Pop III galaxy with nebular 
emission rather than a SMDS (later stellar populations e.g. 
Pop II would also be missing these emission lines, but would 
not be as bright as either Pop HI galaxies or SMDSs). We 
do, however, note one caveat: if there is any supernova (SN) 
explosion that can result from the end of SMDS evolution, 
there might be another way to make He II 1640 radiation. 
When the SN remnant shock reaches the radiatively cool- 
ing stage of its evolution (i.e. when postshock gas cools ra- 
diatively faster than it does by adiabatic expansion), the 
shock becomes a "radiative shock", and that usually means 
that gas cools from a postshock temperature above a mil- 
lion degrees, down to lO^K or below, and He II line emis- 
sion will also occur. The shocks that do this need not only 
be SN explosion shocks, but could also be halo virialization 
shocks, for halos large enough to have virial T high enough 
to ionize He II to He III. On the other hand, as discussed 
above, Heger (personal communication Oct 2011) finds that 
(in the absence of rotation), fusion powered stars more mas- 
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sive than 153, DOOM© collapse directly to supermassive black 
hole seeds rather than blowing up as SN. 

While the detection of emission lines most Hkely in- 
dicates that the object is a Pop III galaxy rather than a 
SMDS, the lack of emission lines leaves both possibilities 
still open. One might therefore ask how well the underly- 
ing continuum spectrum can be determined with JWST. 
The UV continuum slopes for galaxies in the redshift range 
2 — 8 hav e been analyzed using HST data in the liter- 
ature fe.g iBouwens et al] |2009|. i2010l : iDunlop et all l201ll : 



iFinkelstein et al.ll201ll : iMcLure et all 120111 ). The value of 
P can be determine d by converting phot om etric colors in 
releva nt filters (as in lBouwens et aL ll201(t) orlDunlop et al.l 
l|201lf ) for HST), but as noted in lMcLure et al.l (|201lh the 
photometric errors can be quite large, of 0(1). Based on our 
initial estimates it will be difficult to disentangle the SMDS 
from PoplII galaxies based on UV continuum slopes calcu- 
lated from AB magnitude colors in NIRCam or via spec- 
troscopy with NIRSpec. A detailed study of how well this 
separation can be done based on UV spectra is the subject 
of future work. 

The Hell line in Pop HI type A galaxies due to nebular 
emission at 0.1640 iJ.m would fall within the F200W filter of 
JWST for sources at redshifts 9.7 < z < 12.7. The strength 
of the line is modeled in the Yggdrassil code for all PopIII 
galaxies we have considered. Since this line is pronounced 
in Pop HI galaxies but not in DS, one could examine the 
difference in the two magnitudes miso — Jn2oo, which would 
be significantly more negative (i.e bluer miz,o — m^oo colors) 
for DS than for the galaxies. One should be able to see this 
effect for objects a.t z < 12. At higher redshifts, however, 
the Gunn-Peterson cuts off significantly the fluxes in the 
F150W filter, so that it would be impossible to distinguish 
an increase of F200W flux (due to Hell in Pop HI galaxies) 
from a decrease of F150W flux (due to Ly-a absorption). 

Another approach to distinguish between different types 
of o bjects is their loc ation in color/color plots. Previ- 
ouslv llnoue et al.l l|201lh and IZackrisson et al.l l|2011al ) used 
color/color plots to distinguish between different types of 
galaxies: ones with Pop III.l stars vs. ones containing a later 
population of stars. In Figure [16] we try out the possibil- 
ity of differentiating DS vs. galaxies, based on their loca- 
tions in color/color plots. In the left and right panels we 
study Pop HI instantaneous burst galaxies of Type A (max- 
imal nebular emission) and Type C (no nebular emission) 
respectively. We focus here on objects at z=12 as this is the 
most favorable redshift to look for SMDSs. The empty cir- 
cle (cross) symbols correspond to magnitudes for SMDS of 
IO^Mq (1O''M0 ). The solid lines represent the evolution- 
ary tracks of Pop HI galaxies obtained using the Yggdrasil 
model grids, with points marking three different ages (di- 
amonds for IMyr, triangles for 3 Myr, and squares for 10 
Myr). We note that, due to the similar temperatures, SMDS 
formed "with capture" of either 10® M© or lO^M© occupy 
the same spot on the diagrams. 

Pop IIIA galaxies with lifetimes less than 10 Myr will 
exhibit redder colors than SMDS in the mase — 771444 (see 
lower left panel) due to the increased fluxes in the F444W 
filter due to the Balmer emission lines from the galaxies. 
One might therefore hope to distinguish between SMDS and 
Pop HI type A galaxies at z ~ 12 would be to look for red 
colors in 777356 — 777444. Indeed for the case in Figure 1161 this 



technique would work: only the Pop III.l galaxies are bright 
enough to compete with SMDS (see Figure llSp . yet these 
are in a distinctly different location on the color/color plot 
from the SMDS. However, in the figure we have taken the 
stellar mass of the galaxies to be 10® M©, while this num- 
ber could be an order of magnitude higher, which would 
drive the Pop HI lines in the Figure closer to the SMDSs. 
As pointed out before, the error bars in the magnitudes for 
lOo" detections are ~ 0.15, comparable to the differences in 
magnitudes between galaxies and SMDS. In addition, in the 
figure we have taken a specific star formation rate (instan- 
taneous burst). For comparison we have also tried a con- 
stant SFR and found that the results do not change much 
(the colors become slightly bluer). Differentiating between 
SMDS and Pop III galaxies (of uncertain stellar mass and 
metallicity) with such color-color plots will be difficult. 

Differentiating between SMDS and galaxies containing 
Pop HI galaxies is an important issue. Using JWST, the 
best bet is to look for emission lines of He 1640 or Ha. If 
these are found the object is not likely to be a SMDS. On 
the other hand, if these are not found, then differentiating 
between SMDS and Pop III.l galaxies may be difficult with 
JWST. Further studies with other instruments, specifically 
ground based spectrometers, may prove to be helpful. 



7 SUMMARY AND CONCLUSIONS 

The first phase of stellar evolution may have been dark 
stars, powered by dark matter annihilation. These form in- 
side early 10® - 10* halos at z=10-50. Initially DS are 
puffy objects with masses of 1 — IOMq and radii ~10 A.U. 
As long as they are DM powered, their surface temperatures 
(~ lO^K) remain cool enough to allow continued growth 
via accretion of baryons until they become supermassive 
MsMDS ~ 1O®,1O''M0. The requisite DM fuel can be ac- 
quired in two ways: (i) extended adiabatic contraction due 
to DM particles on chaotic or box orbits in triaxial haloes 
and (ii) capture of DM particles via elastic scattering off nu- 
clei in the star. In this paper we have studied the detectablil- 
ity of Supermassive Dark Stars formed via both mechanisms 
with upcoming JWST observations. 

In order to determine their observational characteris- 
tics, we obtained the spectra of SMDSs with the TLUSTY 
stellar atmospheres code (Figure [T]). We used N- Body simu- 
lation s for structure formation at high redshifts (|lliev et al.l 
|2010D to obtain estimates for the numbers of DM halos ca- 
pable of hosting SDMS (Figure [2] and Table [1]). Then we 
used HST observations to set limits on their detectability. 

Both 10® and 10^ Mq SMDS could be seen in HST 
data and wou l d be detected as J-band dropouts. Since 
iBouwens et al.l l|201ll l report only one plausible z~10 object 
in the data, we used the fact that at most one observable DS 
at this redshift can exist to obtain bounds on the possible 
numbers of DS in Eqns. [TT1 and [T2l 

SMDSs are bright enough to be seen in all the wave- 
length bands of the NIRCam on JWST, while detection is 
more difficult in the less sensitive higher wavelength MIRI 
camera. We showed that SMDSs could be seen as J-band, 
H-band, or K-band dropouts, which would identify them as 
Z'^-^IO, 12, and 14 objects respectively. 

The strong point of JWST will be its sensitivity to 
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Figure 16. Signatures of SMDS and instantaneous burst Pop III galaxies in m2rr — masg vs m200 — "1277 (top row) and mase — 711444 vs 
m277 — m356 (bottom row) color diagrams. The left column corresponds to Type A Pop III galaxies (maximal nebular emission) and the 
right column to Type C Pop III galaxies (no nebular emiss i on). Th e solid lines are evolutionary tracks for Pop III galaxies obtained using 
the Yggdrasil model grids introduced in lZackrisson et al. I l l2011al ). The points along the evolutionary tracks single out the evolution at 
three different ages of the galaxies. IO^Mq (IO'^Mq ) SMDS are represented by circle (cross) symbols in the diagrams. For the extended 
AC case we chose a larger size symbol compared to the SMDS formed "with capture" . 



longer wavelengths than HST, corresponding to light from 
higher redshifts where SMDSs may be found. While JWST 
is not particularly better than HST at finding J-band 
dropouts, it will be significantly better at finding SMDS as 
H-band and K-band dropouts. 

We can summarize our predictions for the numbers of 
SMDS seen as H-band dropouts with JWST as; 

Nots = 4.4x 10Vsmd./At(e/150arcmin)2 {Mds = W'^Mq) 

(13) 

Nobs = 2.4x loVsmd. /■At(e/150arcmin)' {Mds = Mq) 

(14) 

where we have scaled the results to 150 arcmin^ survey area, 
which would require multiple surveys by JWST. 

Similarly, our predictions for the numbers of SMDS seen 



as K-band dropouts are: 

Nobs = 7xl0Vsmd./At(e/150arcmin)2 {Mds = 10'^AfQ,AC) 

(15) 

Nobs = 120fsmdsfAt{e/150&rcmmf {Mds = W'^Mq). 

(16) 

IO'^Mq SMDS formed via Capture are not detectable. 

Although these numbers are quite large, as we have em- 
phasized throughout it is quite likely that fsMDS /a* << 1. 
If the DS survives to 2; ~ 10, HST observations bound this 
product. Our final predictions for numbers of SMDS that 
could be detected by JWST are found in Tables [3] and U 

Differentiating between SMDS and galaxies containing 
Pop III galaxies is an important issue. Using JWST, the 
best bet is to look for emission lines of He 1640 or Ha. If 
these are found the object is not likely to be a SMDS. On 
the other hand, if these are not found, then JWST will have 
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trouble differentiating between SMDS and early galaxies. 
Thus further estimates are required using instruments such 
as Giant Magellan Telescope, TMT, LSST, and others. 

As argued by Heger (Heger, private communication) in 
the absence of a dark star phase, the characteristic mass for 
big BHs at birth is 153,000 Mq (i.e., once a fusion powered 
star accretes this much mass it can no longer sustain hy- 
drostatic equilibrium and collapses directly to a BH. With a 
dark star phase, the DS could instead grow to a larger mass 
while DM powered, and then collapse directly to a BH; thus 
in this case the BH could be born with larger masses. Future 
observations of large BHs might thus be able to differenti- 
ate someday between formation mechanisms via dark stars 
or fusion powered stars. 

SMDS mass as a function of halo mass: Although we 
have assumed in this paper that DS grow to the point where 
they consume most of the baryons in the haloes that host 
them, one can examine how the results would change if we 
were to stop the growth at a smaller fraction of the total 
baryonic content. For the case of "maximal bounds" we can 
show that the resulting predictions for JWST remain iden- 
tical. For example, the case we considered in the paper of 
10*' Mq SMDS that grew inside ~ lO'^M© haloes, can be 
compared instead to the case of 10^ Mq SMDS that grew 
inside ~ 10* haloes. For the case of "maximal bounds", 
which assumes that HST bounds at z=10 apply directly to 
SMDS at z=12 (i.e. that the SMDS at z=12 survive aU the 
way to z=10), we find that our results are completely un- 
changed. The number of 10* M© haloes is smaller than the 
number of lO'^M© haloes at both redshifts z=10 (so the HST 
bounds are weaker) and at z=12 (where the JWST obser- 
vations are made). Thus the two effects cancel exactly. One 
can see this cancellation in the following way. The num- 
bers of SMDS observable in either HST or JWST are given 
by the same equation, Eq.((9]). The two factors ^^^^ ^^nd 
fsMDs{z = Zstart) chaugc depending on the hosting halo 
mass, but their product remains the same since it is set by 
HST bounds in Eq.([lOl). Thus the numbers with JWST are 
unchanged regardless of halo size. 

The current decade is a time of great excitement in the 
physics community regarding the possibility of detection of 
the dark matter particle. Three approaches are being pur- 
sued in the hunt for Weakly Interacting Massive Particles: 
direct detection (including DAMA, CDMS, XENON, CO- 
GENT, CRESST, ZEPLIN, TEXONO, COUPP, and many 
others worldwide), indirect detection (including PAMELA, 
FERMI, ICECUBE), and colliders (LHC). Many of these ex- 
periments have indeed found hints of a signal, though confir- 
mation in more than one type of detector of the same particle 
remains a goal. Dark Stars offer a fourth possibility for the 
discovery of WIMPs, or of learning about their properties. 
If WIMPs are indeed discovered, then it is very reasonable 
to expect to find Dark Stars in the sky that are WIMP pow- 
ered. It is even possible that the WIMPs have the property 
that they will be seen first by JWST before any other exper- 
iment. In either case the prospect of finding a new type of 
star in the next premier NASA mission is greatly exciting. 
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